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SUMMARY 

A study  is  reported  which  generated  design  charts  and 
developed  associated  boundaries  for  identifying  departure  and 
uncoordinated  roll-reversal  flight  characteristics  as  a func- 
tion of  three  aerodynamic  parameters.  This  information  should 
be  valuable  for  specification,  design  and  evaluation  purposes. 

The  investigation  utilized  a large  angle, six-degree-of- 
freedom  digital  computer  program  to  simulate  the  motions  of 
a fighter  performing  a severe  air  combat  maneuver  for  differ- 
ent combinations  of  rolling  and  yawing  coefficient  and  lateral 
control  characteristics.  It  was  demonstrated  that  a direction- 
ally  stable  airplane  at  high  alpha  is  virtually  immune  to  de- 
parture. For  a directionally  unstable  configuration,  a high 
dihedral  effect  will  also  prevent  departures.  The  existence 
of  proverse  lateral  control  yaw  characteristics  minimizes  the 
values  of  directional  and  lateral  stability  required  to  avoid 
departure  or  uncoordinated  roll  reversal. 
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High  angle-of-attack  capability  has  been  shown  to  signifi- 
cantly enhance  the  air  combat  maneuvering  (ACM)  effectiveness 
of  fighter  airplanes  (e.g.  see  ref.  1) . Unfortunately,  it  has 
not  been  possible  to  design  for  such  a capability  with  any  de- 
gree of  confidence.  This  is  primarily  attributable  to  the  fact 
that  there  is  no  criterion  which  identifies  the  susceptibility 
of  a specific  design  to  depart  from  controlled  flight  when  per- 
forming rolling  maneuvers.  Consequently,  departures  are  first 
discovered  during  flight  tests  of  the  full-scale  article.  At 
this  point  in  the  airplane  development  cycle,  the  only  recourse 
is  to  incorporate  a device  which  alerts  the  pilot  to  an  approach- 
ing dangerous  flight  regime,  or  to  automatically  limit  the  con- 
trol authorities.  Nevertheless,  under  certain  conditions,  pi- 
lots do  inadvertently  enter  the  flight  regime  where  departures 
and,  consequently,  spins  occur.  Also,  an  alpha  limitation  im- 
posed by  departure  susceptibility  denies  full  use  of  a configur- 
ation's inherent  ACM  capability. 

It  should  be  recognized  that  any  airplane  will  begin  to 
rotate  (i.e.  change  heading)  at  high  alpha  and  enter  a large 
angle  motion  if  a yawing  moment  can  be  generated  and  held  suffi- 
ciently long  through  the  application  of  asymmetric  power  or  dis- 
placement of  an  effective  aerodynamic  control.  The  motion  of 
concern,  however,  is  one  in  which  the  existence  or  direction  of 
a yawing  motion  is  contrary  to  what  had  been  anticipated,  and  the 
motion  occurs  so  quickly  that  the  airplane  is  in  an  incipient 
spin  phase  before  the  pilot  can  take  corrective  action.  The  ob- 
jectives of  the  efforts  reported  herein,  therefore,  are  to  devel- 
op design  charts  and  associated  departure  boundaries  with  which 
to  identify  a particular  design's  degree  of  susceptibility  to 
this  departure  phenomenon.  Such  data  would  be  useful  as  an  eval- 
uation, specification  and  design  tool. 
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SYMBOLS 


wing  span 


q free-stream  dynamic  pressure 

S wing  area 

C wing  mean  aerodynamic  chord 

rolling  moment  coefficient,  rolling  moment  (qSb) 
yawing  moment  coefficient,  yawing  moment/ (qSb) 

C.  lateral  stability  derivative,  3C./3B 

C directional  stability  derivative,  3C  /3B 

nB  n 

C yaw  due  to  lateral  control  derivative,  3C  /36_ 

n*  1 n a 

6a 

I ,1,1  moment  of  inertia  about  X,  Y,  and  Z body  axis, 

X Y Z 

respectively 

Ixz  product  of  inertia,  positive  when  the  principal 

X axis  is  inclined  below  the  X body  axis  at  the 
aircraft  nose 

R^  Reynolds  number 

ot  angle  of  attack,  positive  when  X body  axis  is 

above  projected  relative  wind  vector 

$ sideslip  angle,  positive  when  the  relative  wind 

vector  is  to  the  right  of  the  X-Z  plane. 

<$a  lateral  control  deflection,  positive  to  produce 

left  (negative)  rolling  moment 


slugs 

lb/ft3 


per  deg 
per  deg 
per  deg 


slug-f t: 


slug-f t3 
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TECHNICAL  APPROACH 


Introduction 

The  technical  approach  used  for  this  study  consisted  of: 
o selecting  key  aerodynamic  models, 

o analytically  generating  airplane  time  responses  for  combina- 
tions of  these  models  while  performing  a severe  ACM, 
o tabulating  and  analyzing  the  responses, 

o constructing  departure  susceptibility  design  charts,  and 
o developing  departure  boundaries  from  these  charts. 

The  departure  susceptibility  design  charts  were  based  on 
the  results  of  large  angle,  six-degree-of-f reedom  computations, 
since  departure  from  controlled  flight  and  attendant  incipient 
spin  motions  are  large  amplitude,  coupled  motions  which  reflect 
the  influence  of  gyroscopic  and  kinematic  effects.  Airplane 
departures  are  not,  for  example,  the  result  of  a gust  disturb- 
ance which  develops  into  a large  amplitude  motion  through  the 
means  of  an  unstable  oscillation.  The  employment  of  linearized, 
limited  degree-of-f reedom  equations  of  motion,  consequently, 
would  not  identify  the  departure  phenomenon. 

It  was  believed  that  most  of  the  important  aerodynamic 
parameters  involved  in  the  departure  phenomenon  had  already  been 
identified.  For  instance,  a large  angle,  six-degree-of-f reedom 
analytical  study  (ref.  2)  sponsored  by  NASC  in  1967  determined 
that  adverse  yaw  due  to  lateral  control,  dihedral  effect  and 
directional  stability  were  the  critical  aerodynamic  characteris- 
tics involved  in  promoting  a spin.  Analytical  and  flight  test 
investigations  since  then  have  corroborated  these  findings. 

It  was  assumed  that  control  induced  departures  are  func- 
tions of  the  relative  magnitude  of  these  key  aerodynamic  param- 
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eters  and  their  variation  with  angle  of  attack,  as  well  as  the 
airframe's  inertial  characteristics.  Consequently,  this  study 
endeavored  to  determine  how  the  relative  values  of  these  aero- 
dynamic and  inertial  variables  influence  departure  susceptibil- 
ity. The  ability  to  construct  departure  susceptibility  design 
charts  was  predicated  on  ascertaining  that  a reasonable  number 


of  aerodynamic  (C, 


C and  C 
n n. 


) and  inertial  models  could 


represent  the  spectrum  of  fighter-type  airplanes, 


Aerodynamic  Models 

The  aerodynamic  models  discussed  in  this  section  are  pre- 
sented in  Appendix  A.  All  of  the  aerodynamic  parameters  were 
a function  of  angle  of  attack  and,  depending  upon  the  parameter, 
such  other  variables  as  control  deflection  and/or  sideslip  angle. 
Since  supersonic  flight  was  not  investigated,  Mach  number  effects 
were  not  included  in  the  models. 

The  first  phase  of  the  study  was  the  selection  of  the  C^, 

C and  C models  to  be  investigated.  Models  were  desired 
n nfi 

a 

which  would  encompass  all  possible  fighter  configurations.  To 
accomplish  this,  high  angle-of-attack  data  for  fighter-type  air- 
planes were  sought.  The  airplanes  for  which  data  were  obtained 
are  listed  in  Table  I;  some  of  these  data  are  presented  in  plot- 
ted form  in  Appendix  A.  High  R^,  data  were  desired  so  as  to 
insure  that  the  models  chosen  represented  full-scale  airplane 
characteristics.  The  models  were  constructed  on  the  basis  of 
these  data.  Since  all  of  the  data  fell  into  a broad  band 
throughout  the  angle-of-attack  range,  it  was  possible  to  repre- 
sent the  aerodynamic  characteristics  of  all  the  airplanes  with 
a reasonable  number  of  models.  As  shown  in  Table  I,  low  Rg  data 
was  also  obtained.  It  is  interesting  to  note  that  this  data 
fell  within  the  selected  models  for  angles  of  attack  of  40  deg 
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and  above. 

Rolling  Moment  Coefficient 

The  rolling  moment  coefficient,  C^,  was  modelled  as  a func- 
tion of  a and  8.  The  general  trend,  for  a given  level  of  posi- 
tive sideslip,  was  for  to  be  negative  (stable)  at  an  angle 
of  attack  of  zero  degrees  and  become  increasingly  so  with  increas- 
ing alpha  up  to  approximately  the  stall  region,  after  which  no 
further  increase  in  stability  was  realized;  or,  more  generally, 
a decrease  in  stability  was  experienced  for  increasing  alphas. 

In  no  case  was  found  to  become  positive  (unstable) . 

The  models  chosen  for  the  rolling  moment  coefficient  all 
began  at  the  same  value  at  zero  angle  of  attack  and  reached 
their  maximum  stable  value  at  the  same  alpha,  with  one  exception. 
The  exception  was  Model  A,  which  represented  a slatted  configu- 
ration. For  this  case,  , continued  to  increase  in  a stable 
sense  past  the  stall  alpha  of  the  unslatted  configuration, 
reaching  its  maximum  at  a higher  angle  of  attack,  then  return- 
ing to  the  level  of  the  unslatted  at  even  higher  alphas.  The 
models  selected  encompassed  all  of  the  airplane  rolling  coeffi- 
cient characteristics  and,  in  some  instances,  closely  approx- 
imated actual  configurations  (e.g.  F-14,  F-15,  F-lll) . 

Yawing  Moment  Coefficient 

The  yawing  moment  coefficient,  C , is  modelled  as  a func- 
tion of  angle  of  attack  and  sideslip.  As  one  would  expect,  for 
a given  positive  sideslip  angle,  the  airplane  data  shows  to 
be  positive  (stable)  at  low  alphas  until  approximately  the  stall; 
stability  then  decreases  generally  to  an  unstable  value  at  high 
alphas.  This  loss  of  directional  stability  probably  can  be 
traced  to  adverse  sidewash  due  to  a vortex  impinging  on  the 
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vertical  tail  surface  after  wing  panel  stall.  The  level  of  in- 
stability reached  and  the  angle  of  attack  at  which  the  loss  of 
stability  begins  vary  from  configuration  to  configuration.  The 
YA-10  configuration  deviates  from  the  norm  in  that  it  exhibits 
directional  stability  throughout  the  alpha  range. 

The  yawing  momen t models  were  chosen  to  cover  the  band  of 

data,  from  a constant,  stable  C for  all  alphas,  to  a neutrally 

n 

stable  case  (C  equal  to  zero) , to  the  more  common  situation  which 
starts  stable,  then  becomes  unstable  in  the  post-stall  alpha 
range.  (For  this  latter  situation,  only  one  level  of  stability 
below  20  degrees  angle  of  attack  was  investigated  since  the  ACM 
maneuver  employed  herein  exceeded  this  alpha  range  before  a side- 
slip value  was  generated.)  Several  different  combinations  of 
the  level  of  instability  and  the  angle  of  attack  at  which  loss 
of  stability  begins  were  modelled. 


Yaw  Due  to  Lateral  Control 

Yaw  due  to  lateral  control,  C , values  were  chosen  to 

n6 

a 

represent  typical  characteristics  exhibited  by  fighter  airplanes. 

Three  models  were  chosen:  neutral,  (C  identically  zero) ; 

n6 

a 

proverse  (C  remains  negative  throughout  the  alpha  range) ; and 
“«a 

adverse  (C  initially  negative  and  becoming  highly  positive  in 
n6 

a 

the  post-stall  alpha  range) . 


Aerodynamics  Not  Varied 

The  lift,  drag  and  longitudinal  control  power  characteris- 
tics, as  well  as  the  dynamic  derivatives,  were  typical  of  those 
associated  with  many  current  fighter  airplanes  and  were  not  var- 
ied during  this  study.  The  airplane  was  statically  stable  in 
pitch  throughout  the  angle-of-attack  range  and  not  control  limit- 
ed so  that  the  airplane  could  be  trimmed  to  approximately  40  de- 
grees angle  of  attack.  It  should  be  noted  that  C.  was  also 
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not  varied  although  different  levels  of  C were  investigated. 

6a 

Inertia  Parameter  Models 

The  distribution  of  mass  along  the  body  reference  axes  for 
various  airplane  configurations  and  the  models  chosen  are  pre- 
sented in  figure  1.  As  shown,  mass  is  always  distributed  along 
the  fuselage  axis  for  fighters,  but  to  varying  degrees.  This  dis- 
tribution can  be  covered  with  three  inertia  parameter  models. 

Model  A represents  a configuration  which  has  its  mass  distributed 
only  slightly  more  in  the  fuselage  than  in  the  wings.  Model  C 
represents  a configuration  whose  mass  is  concentrated  heavily  in 
the  fuselage.  Since  the  Model  B configuration  lies  between  the 
two  extremes,  it  was  used  extensively  in  the  aerodynamic  parameter 
study.  The  mass  characteristics  for  this  model  are  presented  in 
Table  II. 

In  order  to  determine  the  influence  of  mass  distribution  on 
an  ACM,  the  multiplicative  factors  presented  in  Table  II  were 
used  to  simulate  Models  A and  C.  These  factors  multiplied  only 
the  contribution  of  the  inertia  coupling  (gyroscopic)  terms  in 
the  equations  of  motion.  In  this  manner,  the  gyroscopic  effect  of 
Model  B could  be  magnified  or  diminished  without  having  to  change 
the  inertia  values,  I , I or  I„,  which  would  have  simultaneously 
incurred  unrealistic  corresponding  changes  in  the  aerodynamic 
acceleration  and  damping  characteristics. 

Flight  Condition  and  Maneuver 

The  airplane  was  trimmed  in  a 60-degree  bank  angle  turn  at 
35,000  feet  and  Mach  0.9  before  a rolling  pull-out  type  maneu- 
ver was  performed.  Control  inputs  for  this  maneuver  were: 

o Full  trailing-edge-up  longitudinal  control  deflection 
initiated  at  time  equal  zero  at  a rate  of  30  deg/sec. 

o Full  lateral  control  deflection  initiated  at  time 

equal  1.5  sec  at  a rate  of  30  deg/sec  in  the  direction 
to  unbank  the  airplane. 

o Rudder  remained  undeflected. 

The  controls  remained  fully  deflected  until  eight  seconds  into 
the  time  history,  at  which  time  both  the  longitudinal  and  lat- 
eral controls  were  returned  to  trim  at  30  deg/sec. 
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This  maneuver  was  chosen  because  it  represented  an  ACM 
from  which  departures  are  often  experienced.  It  is  also  severe: 
controls  are  fully  deflected,  the  airplane  is  pitched  up  rapid- 
ly through  the  stall  and,  being  uncoordinated,  large  sideslip 
angles  are  generated. 

Parameter  Investigation 

A matrix  of  the  aerodynamic  and  inertia  parameter  model 
variations  investigated  is  shown  in  Table  III.  To  evaluate 
these  parameter  variations,  20-second  time  histories  of  the  air- 
frame motions  were  generated  in  response  to  the  previously  de- 
scribed control  inputs.  These  motions  were  computed  using  a 
large  angle,  six-degree-of -freedom  digital  computer  program  at 
the  NADC  computer  facility.  The  program  used  non-linear  tabu- 
lated data  for  the  aerodynamic,  atmospheric  and  control  inputs. 
When  required,  aerodynamic  parameter  tables  could  be  programmed 
as  functions  of  up  to  three  independent  variables  (e.g. , angle 
of  attack,  angle  of  sideslip  and  control  deflection) . Control 
deflection  input  tables  were  programmed  as  a function  of  time. 

The  resultant  time  history  output  (i.e.  the  vehicle  re- 
sponse to  the  control  inputs)  were  plotted  on  a CALCOMP  drum 
plotter.  Up  to  64  output  parameters  of  interest  could  be  plotted 
against  time  for  each  computer  run.  For  the  purpose  of  this 
study,  the  following  sixteen  parameters  were  plotted  because 
they  were  deemed  to  be  of  particular  significance: 

pitch  angle  elevator  deflection 

bank  angle  lateral  control  deflection 

yaw  angle  flight  path  angle 

pitch  rate  altitude  change 

roll  rate  range  position 

yaw  rate  dynamic  pressure 

angle  of  attack  velocity 

sideslip  angle  rotation  rate 
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Prodcedure  for  Developing  Design  Charts 

Each  parameter's  influence  on  the  computed  time  history 
was  determined  by  examining  the  path  of  the  center  of  gravity 
and  the  airplane's  motions  about  the  center  of  gravity,  using 
the  traces  cited  in  the  previous  section.  The  following  time 
history  parameters,  considered  to  be  significant  indicators 
of  these  motions,  were  then  tabulated: 

• peak  yaw  rate  magnitude  and  sign 
o time  to  maximum  yaw  rate 

o time  at  which  yaw  rate  exceeded  +1  rad/sec 

• last  alpha  peak  value  prior  to  lateral  control  removal 

• approximate  second  order  damping  ratio  of  alpha  trace 
o time  for  bank  angle  to  reach  zero  deg  (and  whether 

reversal  occurred) 

• incremental  peak  bank  angle  attained  before  lateral 
control  was  removed 

These  parameters,  in  turn,  were  then  reviewed  to  determine  if 
the  airplane  had  departed  from  controlled  flight.  Subjective- 
ly, the  concept  of  airplane  departure  is  easily  defined;  it 
means  the  pilot  has  lost  positive  control  of  the  airplane. 
Although  an  objective  definition  of  departure  was  not  avail- 
able, four  of  the  above  seven  time-history  parameters  (iden- 
tified by  •)  were  found  to  be  the  most  useful  in  identifying 
departure  susceptibility.  These  parameters  were  plotted  for 
each  model  investigated  and  were  used  as  departure  suscepti- 
bility design  charts.  These  charts  were  used  to  construct 
departure  boundaries. 

Presentation  of  Results 

Figure  2 presents  departure  susceptibility  design  charts 
in  terms  of  the  following  time-history  parameters: 

Figure  2a  - last  alpha  peak  prior  to  lateral  control  re- 

i- 

i 

i 
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moval 

Figure  2b  - approximate  second  order  damping  ratio  of 
alpha  trace 

Figure  2c  - peak  yaw  rate 

Figure  2d  - incremental  peak  bank  angle  reached  at  time 
lateral  control  removed 

These  figures  present  the  value  of  the  time  history  param- 
eter of  interest  for  the  combinations  of  aerodynamic  models 

(C  , C0,  C ) shown  by  the  matrix  of  computer  runs  in  Table 
n 4 n6 

a 

III  for  inertia  parameter  Model  B,  excluding  unstable  Cn  models 

(2)  and  (3).  Each  of  the  figures  contains  three  plots:  one 

each  for  the  unstable,  neutral  and  stable  Cn  models  investigated. 

Each  of  these  plots , in  turn , presents  the  value  of  the  time- 

history  parameter  for  each  combination  of  C and  C models 

n x 
^a 

studied  by  presenting  the  time-history  parameters  as  a function 


of  the  C„  model  for  constant  values  of  C 
H n. 


To  facilitate 


interpretation,  the  models  are  spaced  along  the  abscissa 
according  to  the  value  of  the  C slope  of  the  C models  at 

HQ  l 

40  degrees  angle  of  attack,  with  each  model  indicated  by  sym- 
bols according  to  the  key.  The  40  degree  alpha  C»  value  was 

B 

chosen  because:  (1)  the  models  exhibit  constant  vs  a charac- 
teristics for  a given  level  of  sideslip  above  this  angle  of 
attack  and  (2)  because  the  angle  of  attack  is  either  converg- 
ing to  a trim  level  of  approximately  40  degrees  or  diverg- 
ing to  a higher  alpha  at  the  time  the  controls  are  returned 
to  trim. 
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DISCUSSION 


Design  Charts 

Figures  2a  and  2b  show  how  the  angle-of-attack  behavior 
varies  for  different  combinations  of  the  aerodynamic  models. 

For  a stable  Cn  model,  the  angle-of-attack  time  history  traces 
are  damped  and  the  angle  of  attack  approaches  its  trim  value 
of  approximately  40  degrees  regardless  of  the  C£  and  C models 

*^a 

employed.  For  neutral  or  unstable  yawing  moment  models,  C^  be- 
comes the  most  important  parameter  in  determining  if  the  alpha 
time  history  is  convergent  or  divergent.  If  Cn  is  unstable 
above  the  stall,  as  is  the  case  for  most  airplanes,  a fairly 
high  level  of  dihedral  effect  is  required  to  prevent  the  air- 
plane from  diverging  in  angle  of  attack  to  levels  where  spins 
could  occur.  The  C 


n, 


level  also  influences  the  degree  of 


divergence  and  the  magnitude  of  the  peak  alpha  value,  but  to  a 
lesser  extent  than  the  C„  level. 


I Figures  2c  and  2d  show  that  the  direction  in  which  the 

* airplane  rolls  and  yaws  in  response  to  the  control  inputs 

I also  depends  strongly  upon  the  C^  characteristics.  For  a 

stable  yawing  moment  model,  the  airplane  will  roll  in  the  direc- 
tion commanded  (negative)  accompanied  by  a coordinating  (nega- 
tive) yaw  rate,  whereas  the  opposite  occurs  for  an  unstable 
yawing  moment  model.  These  responses  (both  the  favorable 
and  the  undesirable)  are  significantly  minimized  with  increas- 
ing dihedral  effect.  As  becomes  less  proverse,  a favor- 

but,  in  this  case,  an  undesir- 

the  neutral  C model,  the  j 

^ ! 

general,  do  not  vary  much  i 


l 


able  response  is  also  minimized; 
able  response  is  maximized.  For 
peak  bank  angle  and  yaw  rate,  in 


11 


NADC-76154-30 


with  the  C.  level  (for  C greater  than  -.0015),  but  are 

n X, . 


'B  e 

significantly  influenced  by  the  C 


n. 


model.  The  airplane 


a 


rolls  in  the  direction  commanded  and  experiences  a coordinating 
yaw  rate  for  the  proverse  C , whereas  the  opposite  occurs 


for  adverse  C . For  C equal  to  zero,  the  overall  response 
^a  ^a 

in  roll  and  yaw  averaged  out  to  approximately  zero  by  the 
time  the  lateral  control  was  neutralized. 


These  design  charts  show  that  a stable  yawing  moment 
characteristic  is  the  most  effective  parameter  governing  de- 
parture prevention.  It  also  tends  to  prevent  control  rever- 
sal during  an  uncoordinated  rolling  maneuver  at  high  alphas. 

A designer's  recourse  when  confronted  with  a configuration 
directionally  unstable  over  some  alpha  range  is  to  insure  a 
high  level  of  dihedral  effect.  For  this  combination  of  aero- 
dynamic characteristics,  departures  are  eliminated  at  the  ex- 
pense of  inducing  roll  reversal.  This  constitutes  a trade-off 
between  a safety-of-f light  condition  and  a flight  characteris- 
tic that  might  be  regarded  as  annoying  (i.e.  use  of  rudder 
pedals  required) . It  is  evident  that  these  design  charts  can 
be  used  for  configuration  studies  in  which  directional  stabil- 
ity, dihedral  effect  and  lateral  control  characteristics  may 
be  traded-off  for  the  high  alpha  maneuvering  responses  desired. 

The  design  charts  in  figure  2 present  only  the  models 
which  differed  at  high  angle  of  attack  (>.  40  degrees)  . The 
influence  of  the  rate  of  onset  of  these  high  alpha  values 
are  not  shown.  This  situation  was  investigated  for  the  un- 
stable C and  adverse  C models.  Figure  3 presents  the 

6 

a 

same  four  time-history  parameters  as  those  presented  in  figure 
2,  plotted  against  the  angle  of  attack  at  which  the  model 
assumed  its  high  alpha  value  for  a given  sideslip  angle.  As 
can  be  seen,  the  low  alpha  portion  of  the  curve  does  in- 
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fluence  the  magnitude  of  the  angle-of -attack  characteristics, 
bank  angle  and  yaw  rate,  but  it  does  not  change  the  results: 
angle  of  attack  still  diverges  for  all  cases  and  the  airplane 
experiences  a roll  reversal  for  this  uncoordinated  maneuver. 

A departure  boundary  based  solely  on  figure  2 data  would, 
therefore,  still  be  valid  despite  exclusion  of  this  model 
characteristic. 

Departure  Boundary 

Yaw  rate  alone  is  not  a good  indicator  of  departure,  since 
many  combinations  of  the  aerodynamic  parameters  investigated 
result  in  yaw  rates  that  are  opposite  (unfavorable)  to  the 
commanded  lateral  control  displacement,  but  are  accompanied 
by  a roll  reversal.  Consequently,  the  airplane  yaws  in  the 
same  direction  it  rolls,  and  contrary  to  command.  This  roll 
reversal  is  not  a departure;  it  is  actually  the  "safety  valve" 
that  precludes  departure  and  presents  a strong  signal  to  the 
pilot  to  enter  the  control  loop  and  coordinate  the  maneuver. 

Departures  leading  to  incipient  spins  have  a common 
characteristic:  the  angle  of  attack  must  be  sustained  above 

the  high  alpha  trim  point.  Hence,  the  alpha  time-history 
parameters  were  used  to  determine  if  departure  had  occurred. 

The  model  used  for  the  study  trimmed  at  an  alpha  of  approx- 
imately 40  degrees.  Departure  was  assumed  if  the  alpha  time 
history  demonstrated  progressively  higher  peak  alphas  from 
trim. 


The  angle-of-attack  information  contained  in  figures  2a 
and  2b  were  cross-plotted  as  functions  of  C and  C.  for  each 


n; 


model. 


B 


'6 


C was  the  slope  of  the  C models  at  40  degrees 
nB  n 
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alpha  and  was  used  for  the  same  reasons  previously  discussed 

relative  to  C.  . Values  of  C and  C.  which  resulted  in 

ne 

attaining  given  levels  of  peak  alphas  were  obtained  from  these 
cross  plots  and  were  plotted  as  boundaries,  as  illustrated  in 
figure  4 for  the  adverse  C model.  Figure  4 also  shows  the 

% 

C and  C0  combinations  which  resulted  in  constant  alpha  trace 

ne  p 

damping  ratios.  In  order  to  choose  a departure  boundary  which 
considered  both  a divergent  oscillation  (i.e.,  negative  damping) 
situation  and  an  alpha  build-up  to  high  values,  although  oscilla- 
tions were  damped,  it  was  decided  to  use  the  50  degree  alpha 
peak  (10  degrees  above  trim)  curves  in  conjunction  with  the  zero 
damping  ratio  curves,  choosing  the  most  conservative  as  the 
departure  boundary. 


It  should  be  noted  that  the  aerodynamic  model  employed 
for  this  study  consisted  of  static  aerodynamics  and  dynamic 
derivatives.  For  this  reason,  conclusions  concerning  whether 
or  not  a spin  occurs  should  not  be  drawn  from  these  departure 
boundaries.  To  predict  spins,  it  is  necessary  to  model  rotary 
aerodynamics,  as  has  been  shown  in  ref.  3. 


A departure  boundary  for  each  level  of  C is  shown  in 

n6 

a 

figure  5.  For  combinations  of  C and  C above  the  boundary, 

ne 

departure  does  not  occur.  For  combinations  below  the  boundary, 
departure  occurs  for  the  maneuver  used  in  this  study.  For  less 

severe  maneuvers,  departure  may  not  occur  even  for  C and 

nB 

C.  combinations  below  the  boundary,  but  the  danger  of  depar- 

ts 

ture  would  be  present.  As  an  example,  several  computer  simu- 
lations were  made  for  which  the  controls  were  removed  at  times 
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less  than  eight  seconds  (e.g.  two  and  four  seconds) . For  very 
short  control  input  durations,  departure  would  not  occur  for 
any  combination  of  aerodynamic  models  investigated,  since  the 
airplane  was  unable  to  develop  any  sizable  motions.  However, 
as  the  input  duration  increased,  departures  would  be  experi- 
enced for  the  aerodynamic  parameter  combinations  indicated  by 
the  boundaries. 


The  departure  boundaries  show  that  for  large  stable  values 
of  C.  , fairly  large  unstable  C values  can  be  tolerated  with- 

£e  ne 

out  departure  susceptibility,  regardless  of  C characteristics. 

n6 


a 

For  low  levels  of  C.  , more  stable  (or  at  least  less  unstable) 

8 

C values  would  be  required  for  a configuration  to  remain  depar- 
nB 

ture  resistant.  For  an  airplane  which  exhibits  a positive 

(stable)  C , the  influence  of  C.  becomes  unimportant  except 
nB  *8 

for  a configuration  with  virtually  no  dihedral  effect. 


1 
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The  results  presented  in  figures  2 through  5 were  derived 

using  the  inertia-parameter  Model  B.  A limited  number  of  time 

histories  were  generated  for  the  inertia-parameter  Models  A and 

C.  In  figure  6,  the  angle-of-attack  design  chart  parameters 

obtained  for  inertia-parameter  Model  B are  compared  with  the 

values  obtained  with  Models  A and  C.  These  data  are  presented 

since  they  were  employed  to  develop  the  departure  boundaries. 

As  shown,  significant  differences  were  obtained  only  for  Model 

C at  some  C values.  Even  these  differences,  however,  would 
p 

not  significantly  affect  the  departure  boundaries  which  were 
based  on  zero  damping  or  a 50  degree  alpha  peak  value.  Ad- 
mittedly, the  inertia-parameter  model  study  was  limited  and 
further  study  would  be  warranted  to  fully  document  the  impor- 
tance of  mass  distribution  on  the  departure  boundaries. 
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The  effect  of  maintaining  directional  stability  to  a higher 
alpha,  unstable  model  (2)  , is  shown  in  figure  7.  The  de- 
sign chart  parameters  show  that  maintaining  directional  sta- 
bility to  a higher  angle  of  attack  results  in  lowered  peak 
alphas  and  an  increased  damping  ratio,  with  an  increased  peak 
yaw  rate  and  bank  angle.  As  shown  in  figure  8,  a sizable 
shift  in  the  alpha  at  which  directional  stability  is  lost  does 
not  change  the  slope  of  the  boundary  and  results  in  only  a 
small  shift. 

Uncoordinated  Roll  Reversal  Boundary 


Examination  of  the  design  charts  (figure  2)  showed  that 

for  a large  group  of  the  cases  studied  the  airplane  rolled 

in  a direction  opposite  to  that  commanded.  The  initial  motion 

was  in  the  commanded  direction,  to  varying  degrees;  it  then 

reversed  and  rolled  in  the  opposite  direction.  This  maneuver 

is  uncoordinated  and,  therefore,  large  sideslip  angles  are 

generated  at  large  alpha.  A commanded  negative  roll  rate  at 

positive  angles  of  attack  produces  a sideslip  angle  which, 

for  an  unstable  C , will  introduce  a positive  yaw  rate.  This, 

nB 

in  turn,  leads  to  more  negative  sideslip  and  a positive  roll- 
ing moment,  due  to  dihedral  effect  (negative  C ) . An  adverse 

p 

C aggravates  the  situation  by  producing  additional  positive 
"6a 

yaw  rate  due  to  the  deflection  of  the  lateral  control.  For 
a stable  C , this  effect  is  combatted  since  the  magnitude 

ne 

of  the  negative  sideslip  is  depressed  or  positive  sideslips 
are  realized. 

It  is  possible  to  generate  boundaries  using  the  bank  angle 
information  in  figure  2 in  the  same  manner  as  was  done  for  the 
departure  boundaries.  These  are  shown  in  figure  9.  It  is  em- 
phasized that  these  boundaries  are  not  departure  boundaries 
but  are  uncoordinated  roll-reversal  boundaries  requiring  the 
pilot  to  apply  rudder  to  coordinate  the  maneuver  or  to  depress 
the  angle  of  attack. 
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CONCLUSIONS 


Departure  susceptibility  and  uncoordinated  roll-reversal 
boundaries  have  been  developed  from  design  charts  generated  by 
considering  a series  of  roll  and  yaw  coefficient  characteris- 
tics for  different  types  of  yaw  due  to  lateral  control  deflec- 
tion. For  ease  of  use,  these  boundaries  are  presented  in  terms 


of  derivatives,  i.e.  C.  , C 

*8  n 


B 


n , 


The  boundaries  and  design 


charts  can  serve  as  valuable  evaluation,  specification  and  de- 
sign tools.  They  were  used  to  provide  the  following  conclu- 
sions : 

o positive  directional  stability  at  high  alpha  prevents 
the  airplane  performing  a severe  ACM  from  diverging  to 
alphas  above  trim  and  into  the  region  where  spins  are 
encountered,  regardless  of  the  level  of  dihedral  effect 
and  the  value  of  C 


n 


6 


o For  directionally  unstable  configurations,  a large 

dihedral  effect  (negative  C.  ) can  avert  departure. 

8 

However,  this  may  require  rolls  to  be  coordinated  to 
avoid  roll  reversal. 

o The  presence  of  a proverse  Cn  characteristic  mini- 
mizes the  value  of  C and  C.  required  to  avoid  de- 

n8  *8 

parture  or  uncoordinated  roll-reversal. 


As  in  all  studies,  the  scope  of  the  effort  was  influenced 
by  the  time  and  resources  available.  The  limitations  placed 
on  the  present  study,  however,  did  not  compromise  the  applica- 
bility of  the  results  and  the  development  of  boundaries  and 
definitive  conclusions.  Nevertheless,  the  following  limita- 
tions should  be  noted.  Variations  of  the  mass-inertial  param- 
eters within  the  range  that  characterize  fighter  airplanes 
seemingly  do  not  change  the  boundaries  which  were  determined. 
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Further  study  is  required  to  make  this  conclusion  unequivo- 
cally. The  aerodynamic  parameters  C and  C , while  not 

m m0 

q p 

treated  herein,  may  cause  displacement  of  the  boundaries. 
Also,  the  boundaries  were  developed  for  airplanes  that  are 
stable  in  pitch  throughout  the  angle-of-attack  range  (between 
0 and  90  degrees)  and  have  no  limitation  placed  on  the  avail- 
able longitudinal  control  authority.  Pronounced  changes  in 
these  characteristics  would  patently  shift  these  boundaries. 
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TABLE  I.-  AIRPLANES  FOR  WHICH  DATA  WAS  AVAILABLE 


High  Re  data  - 13  configurations 


F-4E 

(slatted 

and  unslatted) 

F-5E 

F-14 

(slatted 

and  unslatted) 

YA-10 

F-15 

(Republic 

proposal) 

TA-4 

F-15 

F-18 

F-lll 

T-2C 

F-5A 
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TABLE  II.-  MASS  CHARACTERISTICS  OP  INERTIAL  MODELS 

INERTIAL  MODEL  B 

Ix=  25,000  (Ix-Iy)/mb2=  -0.067  j 

j 

Iy=  135,000  <Iy-Iz)/mb2=  -0.012  ! 

Iz=  155,000  (Iz-Ix)/mb2=  0.079 

xxz-  0 

m=  1025 

INERTIAL  MODEL  A 

(IX“IY)  = °-407  (IX"IY) 

A B 

(IY’IZ)  = 1'°  <1y“IZ) 

A B 

(W  = °-438  'W 

A B 


INERTIAL  MODEL  C 
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(a)  Last  alpha  peak  prior  to  lateral  control  removal  as  a function  of  the  aerodynaaic  models. 
Figure  2.-  Design  charts  for  determining  departure  susceptibility. 
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Incremental  peak  bank  angle  reached  at  time  lateral  control  removed 
Figure  2.-  Concluded. 
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Figure  4.-  Constant  alpha  peak  and  damping  ratio  curves  derived 
from  design  charts  as  a function  of  lateral  and  directional 
stability  for  adverse  C 
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APPENDIX  A 

The  C0 , C and  C models  investigated  in  this  study,  as 
£ n n ft 

6a 

well  as  a compilation  of  the  airplane  data  upon  which  these 
models  were  based,  are  presented  in  this  appendix.  Also  shown 
are  the  aerodynamic  derivatives  which  were  held  constant  during 
the  investigation. 

C„  and  C values  obtained  during  high  Reynolds  number  wind 

Ar  11 

tunnel  tests  of  various  fighter  type  configurations  are  pre- 
sented in  figures  Al  and  A2,  respectively,  as  a function  of  angle 
of  attack  for  10  degrees  of  sideslip.  These  are  sample  plots  of 
the  collected  data  employed  to  construct  the  and  Cn  models. 

Figures  A3  and  A4  show  the  selected  models  as  a function 
of  angle  of  attack  for  10  and  25  degrees  of  sideslip,  respective- 
ly. is  zero  throughout  the  angle-of -attack  range  for  zero 

sideslip  and  is  assumed  to  remain  constant  with  sideslip  for 
sideslip  angles  greater  than  25  degrees. 

C models  are  shown  in  figures  A5  and  A6  for  10  and  25  de- 
grees of  sideslip,  respectively,  as  a function  of  angle  of 
attack.  As  for  the  models,  is  zero  for  zero  sideslip,  and 
the  level  is  constant  for  sideslip  angles  greater  than  25  de- 
grees. 

The  C models  are  shown  in  figure  A7  to  be  a function  of 
n, 

a 

angle  of  attack  and  represent  typical  levels  of  proverse  and  ad- 
verse yaw  realized  with  fighter  airplane  configurations. 

I 

t 

The  aerodynamic  parameters  which  were  not  varied  during  j 

the  study  are  presented  in  figures  A8  through  A24.  As  shown,  * 


A-l 
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normal  and  axial  coefficients,  figures  A8  through  All,  were 
programmed  as  functions  of  angle  of  attack,  sideslip  angle  and 
longitudinal  control  deflection.  The  pitching  moment  coefficient, 
figures  A12  and  A13,  was  also  programmed  as  a function  of  angle 
of  attack,  sideslip  angle  and  control  deflection.  The  pitching 
moment  is  typical  of  what  might  be  found  on  those  fighter  air- 
planes which  are  statically  stable  in  pitch  throughout  the  angle- 
of -attack  range.  The  side  force  coefficient,  C^,  also  programmed 
as  a function  of  a,  8 and  longitudinal  control  deflection,  is 
shown  in  figures  A14  and  A15.  The  lateral  control  derivatives 
were  programmed  as  functions  of  angle  of  attack  and  lateral  and 
longitudinal  control  deflections.  The  longitudinal  control  de- 
flection dependence  is  due  to  the  fact  that  many  current  fighter 
configurations  utilize  differential  elevator  deflection  for  lateral 
control.  The  C.  model  is  shown  in  figures  A16  and  A17,  and  the 

C model  is  shown  in  figures  A18  and  A19.  The  dynamic  deriv- 


atives, C , C , etc.,  are  shown  in  figures  A20  through  A24; 
n m 
P <3 

each  is  a function  of  angle  of  attack  and  represents  characteristic 
values  for  fighter  type  airplanes. 
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Figure  A2.- Yawing  moment  coefficient  as  a function  of  angle  of 
attack  for  various  airplanes.  6-  10  deg. 
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Figure  A3.  - Rolling  moment  coefficient  models  for  0=  10  deg. 


Rolling  moment  coefficient  models  for  3=  25  deg 


Figure  A8.- Normal  force  coefficient  as  a function  of  angle  of  attack 
end  sideslip  angle.  i = 0 deg. 
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Figure  A12 Pitching  moment  coefficient  as  a function  of  angle  of 
attack  and  sideslip  angle,  i - 0 deg. 
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Figure  A16.-  Rolling  moment  due  to  lateral  control  derivative  as  a function 
of  angle  of  attack  and  control  deflection,  i = 0 deg. 
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force  due  to  lateral  control  derivative  as  a function  of 
le  of  attack  and  control  deflection,  i = -30  deg. 


Figure  A21.-  Damping  derivative,  C.  , as  a function  of  angle  of  attack. 


Figure  A22.-  Cross  derivative,  C , as  a function  of  angle  of  attack 


Figure  A23.-  Cross  derivative,  Cf  , as  a function  of  angle  of  attack 


